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ABSTRACT 

Objectives: Respiratory diseases are ranked second 
in Europe in terms of mortality, prevalence and costs. 
Studies have shown that extreme heat has a large 
impact on mortality and morbidity, with a large relative 
increase for respiratory diseases. Expected increases in 
mean temperature and the number of extreme heat 
events over the coming decades due to climate change 
raise questions about the possible health impacts. We 
assess the number of heat-related respiratory hospital 
admissions in a future with a different climate. 
Design: A Europe-wide health impact assessment. 
Setting: An assessment for each of the EU27 
countries. 

Methods: Heat-related hospital admissions under a 
changing climate are projected using multicity 
epidemiological exposure-response relationships 
applied to gridded population data and country-specific 
baseline respiratory hospital admission rates. Times- 
series of temperatures are simulated with a regional 
climate model based on four global climate models, 
under two greenhouse gas emission scenarios. 
Results: Between a reference period (1981-2010) and 
a future period (2021-2050), the total number of 
respiratory hospital admissions attributed to heat is 
projected to be larger in southern Europe, with three 
times more heat attributed respiratory hospital 
admissions in the future period. The smallest change 
was estimated in Eastern Europe with about a twofold 
increase. For all of Europe, the number of heat-related 
respiratory hospital admissions is projected to be 
26 000 annually in the future period compared with 
11 000 in the reference period. 
Conclusions: The results suggest that the projected 
effects of climate change on temperature and the 
number of extreme heat events could substantially 
influence respiratory morbidity across Europe. 
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INTRODUCTION 

Respiratory diseases are ranked second in 
Europe in terms of mortality, prevalence and 
costs. 1 This burden is expected to increase, 
partly due to a changing climate. 2 An environ- 
mental factor with a large impact on mortality 



ARTICLE SUMMARY 



Article focus 

■ To assess how heat-related respiratory hospital 
admissions in Europe could change with increas- 
ing temperatures in the near future. Within a 
range of climate change projections, explore dif- 
ferent impact estimates. 

Key messages 

■ The 30-year mean annual increase in heat-related 
respiratory hospital admissions between the 
periods 1981-2010 and 2021-2050 can be 
counted in the 10s of thousands because of pro- 
jected temperature increases. The increase will 
be different in different parts of Europe with the 
largest relative increase in southern Europe. 

Strengths and limitations of this study 

■ The study takes the spatial variation in climate and 
population density into account in impact calcula- 
tions, making the estimates valid for countries with 
large national differences in these factors. 

■ The study results exemplify some of the variation 
of health impacts under different climate change 
scenarios, as well as the variations within one 
scenario, with different underlying climate 
models. 

■ The use of spatial population data limited the 
study to exploring the change in all age respira- 
tory hospital for studying heat-related illness; 
combined with the future increase in the propor- 
tion of elderly people, this suggests that the 
estimated impacts are likely to be conservative. 



and morbidity in Europe is extreme heat, with 
large effects on respiratory diseases. 3-5 The 
physiological effects of exposure to heat can be 
directly heat-related (heat stroke, heat fatigue 
and dehydration) or can contribute to a wor- 
sening of respiratory and cardiovascular dis- 
eases, electrolyte disorders and kidney 
problems. 4-7 The reasons for an increase in 
respiratory admissions may be several. The 
elderly with respiratory diseases such as 
chronic obstructive pulmonary disease 
(COPD) are less fit and suffer often from 
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circulatory problems. Heat influenced admissions due to 
chronic airway obstruction and asthma increased more 
than admissions due to chronic bronchitis in a study from 
New York City, but the daily number of admissions seldom 
allow a study of specific diagnoses. In a recent expert elicit- 
ation among European researchers engaged in environ- 
mental medicine or respiratory health, extreme heat stood 
out as the most important climate-related pathway to 
adverse impacts on respiratory health, more important 
than changes in air pollutants and allergens. 9 A review 
found that heatwaves have a stronger relative impact on 
mortality than on emergency room visits and hospital 
admissions, suggesting that many individuals die before 
they can get to the hospital. 10 However, several studies 
confirm that heat also affects healthcare utilisation. 3 11 
During the heatwave of 2003 in France, the Assistance 
Public Hopitaux de Paris recorded 2400 additional visits to 
the emergency care units and 1900 additional hospital 
admissions. 12 During the heatwave in California in 2006, 
there were almost 1200 excess hospital admissions together 
with more than 16 000 excess visits to the emergency 
departments. 

It is possible that hospital admissions could increase in 
the future owing to temperature changes projected with 
climate change, assuming no additional acclimatisation. 
Projected changes include an increase in the global 
mean temperature of 1.8 to 4°C by the end of this 
century, with larger increases over land areas at high lati- 
tudes, changing seasonal temperatures and increases in 
the frequency, intensity, duration and spatial extent of 
heatwaves. 13 Recent projections for 2070-2100 suggest 
that maximum temperatures experienced once every 
20 years during the period 1961-1990 could be expected 
as often as every year in southern Europe and every 
third and fifth year in northern Europe. 14 This means 
that Europe should expect an increase in both mean 
temperatures and the number of extreme heat events. 
Recent heatwaves are consistent with projections. For 
example, the heatwave that occurred in Europe 2003 
could be expected to return every 46 000 years, based 
on the temperature distribution for the years 1864- 
2000. 15 However, such extreme hot conditions may 
become much more common at the end of this century 
due to anthropogenic climate change. 16 The uncertainty 
in this estimated return time is quite high, with a lower 
bound of the 90% CI of 9000 years. Even so, another 
heatwave occurred in 2006 with the most anomalous July 
temperature ever measured in Europe. 17 In 2010, 
Eastern Europe experienced a heatwave with summer 
temperatures higher than in the last 140 years, resulting 
in roughly 55 000 excess deaths in Russia. 18 

The possible impacts of these changes on morbidity, 
particularly respiratory health, are relatively unexplored. 

AIM 

The aim of this study is to assess the extent to which 
changes in the frequency of hot days due to climate 



change over the next 40 years could affect heat-related 
respiratory hospital admissions (RHAs) in Europe. Using 
a range of climate projections, we estimate the change 
in hospital admissions between a reference period and a 
future period. 

METHODS 

Climate change and temperature modelling 

The Rossby Centre regional atmospheric climate model 
RCA3 19 was developed by SMHI (the Swedish meteoro- 
logical and hydrological institute) to dynamically down- 
scale results from the global climate models CCSM3, 20 
ECHAM5, 21 HadleyCM3 22 and ECHAM4 23 to a higher 
resolution over Europe. RCA3 is run on a horizontal 
grid spacing of 0.44° (corresponding to approximately 
50 km) and a time step of 30 min. Projections are based 
on the global greenhouse gas emission scenarios A1B 
and A2 (scenarios are described in detail in the 
SRES 24 ). Both have been used in climate change health 
impact assessments (HIA). 25-28 A1B is a 'middle of the 
road' scenario and A2 is considered a high-emission 
scenario, although recent greenhouse gas emissions 
have been higher. Data from one climate model, under 
one climate change scenario, are referred to as one 
climate change projection. 

We use aggregated daily projections of maximum tem- 
perature and relative humidity data to estimate exposure 
for the periods 1981-2020 (reference period) and 
2021-2050 (future period). 

Population and morbidity rates 

The annual country-specific rate of RHAs (ICD 
(International classification of diseases)-9 : 460-519) 
between 2005 and 2010 were extracted from the 
WHO's European Health for All Database (http:// 
data.euro.who.int/hfadb) for the EU27 countries. 
In the dataset, national hospital admissions data are 
provided by the national public health institutes, 
health ministries or corresponding functions. The 
ICD-9 codes were preferred because the epidemio- 
logical studies in the PHEWE (Assessment and 
Prevention of Acute Health Effects of Weather 
Conditions in Europe) study were based on ICD-9 
RHAs. The mean value over the 6 years was used as a 
baseline morbidity rate. In order to have a fine spatial 
resolution, official population data were from the 
History Database of the Global Environment (HYDE) 
theme within the Netherlands Environmental 
Assessment Agency. 29 These data are gridded on a 
0.0833° resolution (approximately 9.45 km) and 
matched with the climatic data by summing the popu- 
lation within each climatic grid cell. 

Exposure-response assumption 

The impact calculations were based on the relationship 
between heat and RHAs (ICD-9 : 460-519) estimated 
within the European PHEWE project. 3 This relationship 
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was based on the analyses of 12 European cities over 
the 1990s, using the daily maximum apparent tempera- 
ture (AT), where AT is a combination of the measured 
temperature and the dew point temperature: 

AT = -2.653 + 0.944T + 0.0153(DT) 2 

where T is the air temperature and DT the dew point 
temperature, which can be derived from temperature 
and relative humidity. 30 AT accounts for heat stress 
related not only to the absolute temperature but also to 
the saturation of the surrounding air, which makes it 
harder to regulate body temperature by sweating, the 
most important mechanism to maintain a healthy body 
temperature. 31 32 The PHEWE model also included 
potential confounders such as air pollution, holidays, 
weekdays, etc. The study used a 0-3 day lag of the 
maximum AT and concluded that the 90th percentile of 
this exposure variable for the summer months, April to 
September, is the appropriate threshold value for the 
heat-morbidity function. 3 We calculated the 90th per- 
centile for the summer months for each grid cell from 
the climate data for the 1990s, yielding individual thresh- 
olds for each grid cell for each climate projection. 

The PHEWE study calculated the relative risk (RR) coef- 
ficients for each of the 12 cities and then combined these 
into two metacoefficients, for north-continental cities 1.012 
(95% CI 1.001 to 1.022) and Mediterranean cities 1.021 
(95% CI 1.006 to 1.036), associated with a 1° increase 
above the temperature threshold. Metacoefficients were 
more suitable as they allowed calculation of the change in 
RHAs for countries with and without a city in the PHEWE 
study. We assigned the Mediterranean coefficient to 
Portugal and the Mediterranean countries. The rest of 
Europe was assigned the RR for the north-continental 
cities. Statistically significant coefficients were found for all 
ages and ages 75+ in the PHEWE study. Because the popu- 
lation data were not stratified by age, we used the PHEWE 
study coefficients for all ages. 

Projections 

The final data used in the impact calculations were 
based on the climate data grid with a resolution of 
0.44°, resulting in 8075 grid cells over Europe. The 
population data, RR coefficients and baseline morbidity 
rates were projected for each grid cell. The countries 
within EU27 were grouped into four regions, northern, 
western, eastern and southern Europe, according to the 
United Nations classification scheme (http://unstats.un. 
org/unsd/methods/m49/m49regin.htm#europe). 

Thus, each grid cell includes (1) a daily time-series of 
0-3 day lag maximum AT for each climate change pro- 
jection, (2) a grid cell-specific temperature-morbidity 
threshold and a location-based RR coefficient and 
(3) population. For each grid cell, the AT was compared 
to the specific threshold for each day in the two time 
periods and risk estimates were calculated for each day. 



The expected number of daily RHAs in each cell was cal- 
culated using the population and the expected daily 
number of RHAs per capita, based on the national 
average of the grid cell. These were combined to calcu- 
late the excess number of RHAs as 

RHA; t = RR ( AT V- Thres ') x Popi x RHApC; 

where RHAi t is the number of RHAs attributable to heat 
in grid cell i at time t, RRj the RR coefficient in grid 
cell i, ( AT; t — ThreSi) + the difference between the 0-3 lag 
AT and the threshold for grid cell i at time t if ATi t were 
greater than the threshold and 0 otherwise, Pop; and 
RHApCi the population and expected number of RHAs 
per capita in grid cell i. 

For each country, the estimated number of RHAs 
attributed to heat was calculated for each climate 
change projection for both the reference period and 
the future period. The estimated number of RHAs 
attributed to heat was then transformed into the propor- 
tion of the expected annual number RHAs for each 
country (data available from the authors). Mean esti- 
mates were calculated for each region for each climate 
change projection. 

Sensitivity analysis 

Because the study results could heavily depend on the 
RR coefficients and thresholds used, we assess whether 
the estimates appeared sensitive to the region-specific 
metacoefficients. We investigated the grids that con- 
tained the cities in the PHEWE study using both the 
location (Mediterranean or north-continental) and city- 
specific RR coefficients. 3 We then compared the results 
with respect to the attributed proportion of RHAs and 
the proportional change in heat-related RHAs. 

RESULTS 

The periods of warm days will increase in the future. For 
the cities included in the sensitivity analysis, the tem- 
perature threshold was exceeded annually, on average, 
by 20 days in the baseline period and 40 days in the 
future period. 

In the future period, approximately 0.4% of the 
annual numbers of RHAs in Europe was estimated to be 
due to heat (table 1), based on the mean estimates over 
the climate change projections. In absolute terms, 
assuming all else as equal, this represents about 26 000 
cases annually in Europe. This should be compared to 
the reference period where approximately 0.18% of all 
RHAs were attributed to heat or about 11 000 cases 
annually. Thus, the results suggest more than a relative 
doubling of the RHAs attributed to heat in Europe. 

On the regional level, the five projections estimate 
increases in the number of heat-related RHAs for 
Europe (table 2). However, in one climate change pro- 
jection, the Czech Republic, Hungary Poland and 
Slovakia were estimated to have a decrease in 
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Table 1 The estimated proportion of RHAs attributed to heat for each region. Intervals describe the highest and lowest 
national estimates in each region 



Region 


1981-2010 


2021-2050 


Change 


Eastern Europe 


0.17% (0.16-0.19%) 


0.31% (0.29-0.35%) 


0.14% (0.11-0.17%) 


Northern Europe 


0.13% (0.10-0.15%) 


0.27% (0.19-0.32%) 


0.14% (0.07-0.17%) 


Southern Europe 


0.23% (0.18-0.26%) 


0.64% (0.42-0.68%) 


0.41% (0.23-0.44%) 


Western Europe 


0.18% (0.16-0.20%) 


0.39% (0.34-0.45%) 


0.21% (0.17-0.26%) 


EU27 


0.18% (0.10-0.26%) 


0.40% (0.19-0.68%) 


0.21% (0.07-0.44%) 



heat-related RHAs (figure 1). The countries with the 
highest estimated increase also show the largest range 
between the highest and lowest estimates. The 
Scandinavian and Baltic countries show the smallest 
range between the highest and lowest estimates along 
with small increases of the mean estimates. 

There is variation among countries, with the largest 
increases in the southern European countries and the 
smallest increase in Eastern Europe (table 2). The rela- 
tive change in the burden of RHAs in relation to the 
climate change scenarios investigated indicates a larger 
relative increase in Mediterranean countries (approxi- 
mately three times) compared to northern European 
countries (approximately two times). 

Sensitivity analysis 

The calculations using the city-specific RR coefficients 
yielded a different proportion and number of estimated 
RHAs than the ones using the two metacoefficients, as 
expected. However, the relative changes in the number 
of hospital admissions attributed to heat between the 
two time periods were not affected by the change in the 
RR coefficients. 



DISCUSSION 

The analysis estimates that the number of RHAs will 
increase in a warmer climate with more hot days. 
Heat-related RHAs were projected to increase two to 
three-fold in the future owing to climate change. 
However, the proportion of respiratory admissions attrib- 
uted to heat would remain rather small. This projection 
is in line with the results of a recent study that estimates 
that respiratory admissions due to excessive heat in 



New York State will increase 2-6 times from the period 
1991-2004 to the period 2080-2099. 33 

As most heat-related health outcomes occur during the 
warmest period of the year, presenting the increase as a 
change in the proportion of the total number of annual 
RHAs can underestimate the additional burden on the 
healthcare system during summer. As the threshold was 
exceeded for approximately 40 days in 1 year, an annual 
increase of 0.21% would result in a 1.9% increase on 
average during these 40 days. The annual numbers are 
used because the available baseline rates of RHAs are 
expressed as an annual average. Applying results to subna- 
tional scales could be inappropriate because national 
averages summarise over considerable heterogeneity. 

The results suggest a larger impact from heat in south- 
ern Europe in the future period, centred on year 2035, 
than in the eastern and northern parts. This is in line 
with many climate change projections showing a larger 
relative increase in the number of extremely hot days in 
southern Europe compared with northern Europe. 14 
However, to some extent, this might also be explained 
by model bias for northern Europe introduced by the 
RCA3 model, where the model appears to underestimate 
temperature for the warmest days. 19 This temperature 
bias is present in the reference period, future period 
and threshold values; therefore, the estimated numbers 
of RHAs from the two periods within each scenario are 
comparable, but comparisons within the same time 
period across scenarios may be inaccurate. 

As the estimates are based on the population size in 
each grid cell, the added burden will be larger in coun- 
tries with larger increases in temperatures in densely 
populated areas. Because the population of a city is con- 
sidered to be a good predictor of the size of the urban 
heat island, 34 an increase in the population living in 



Table 2 Future increase in heat-related RHAs based on the four climate models, under two emission scenarios, as the 


percentage of the annual expected number of RHAs in each region 








Climate model CCSM3 (%) 


ECHAM5 (%) 


HadCM3 (%) 


ECHAM4 (%) 


ECHAM5 (%) 


Greenhouse gas scenario A1 B 


A1B 


A1B 


A2 


A2 


Eastern Europe 0.32 


0.08 


0.18 


0.12 


0.01 


Northern Europe 0.17 


0.09 


0.14 


0.20 


0.08 


Southern Europe 0.51 


0.29 


0.64 


0.45 


0.14 


Western Europe 0.30 


0.11 


0.26 


0.29 


0.06 


EU27 0.32 


0.13 


0.29 


0.26 


0.07 
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Figure 1 The range of the 
absolute increase in RHAs 
attributed to heat between the two 
periods (1981-2010, 2021-2050) 
as a proportion of the annual 
expected number of RHAs for 
each of the 27 countries. The 
points show the highest and 
lowest estimates from four climate 
models under two emission 
scenarios. 
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urban areas will increase the numbers exposed and the 
temperature to which they are exposed. Heat islands 
increase temperatures in urbanised areas compared to 
surrounding areas, also reducing cooling during the 
night-time. 35 These factors combined are likely to 
magnify the health burden during a heatwave. 36 The 
spatial scale of the climate models makes them unable 
to take the urban heat island effect into account. 
Together with the urbanisation of Europe, this could 
potentially result in the underestimation of actual conse- 
quences/RHAs in the future period, because the same 
population size and composition were assumed in both 
periods. In addition, climate change is likely to increase 
ozone concentrations that would add an additional 
health burden for people at risk of respiratory diseases. 
The number of deaths and RHAs, due to a change in 
ozone, is expected to change in the future. 37 A sensitiv- 
ity analysis in the PHEWE study, however, showed that 
the exposure-response relationship for heat did not 
substantially differ between models taking ozone into 
account and models adjusting for NO2 alone. 3 

We used the same thresholds for the reference period 
and future period to isolate the effect of climate change. 
This would tend to overestimate the increased impact 
because there will undoubtedly be biological and/or 
social adaptations that will reduce future health 
temperature-related burdens. However, given the uncer- 
tainty of such adaptation effects, we choose to not 
incorporate such effects in the projections. Recent 
studies in the USA indicate that heat-related health 
burdens have decreased, 38 indicating that some adapta- 
tion is taking place. An opposite trend appears to be 
occurring in Stockholm, with an increase in the risk of 
mortality associated with heat during the 1990s. 39 A 



study looking at the impact of heat on mortality, before 
and after the implementation of a heat-warning system 
in Italy, shows that the effects of extreme heat can be 
reduced, whereas the effects of moderately increased 
average temperatures remain similar. 40 The magnitude 
and the extent of future adaptation is, of course, highly 
uncertain and will vary between and within countries. 
Nevertheless, cities with higher thresholds seem to have 
higher risk ratio coefficients. 3 41 In an effort to estimate 
future and presumably higher thresholds, one must also 
adjust the risk ratio coefficients. This would result in 
fewer days of elevated risk, but the risk increase on each 
occasion could be higher due to the higher RRs. 

Table 2 shows how the results vary by global climate 
models and greenhouse gas scenarios. This study exem- 
plifies the magnitude of the difference between projec- 
tions made by a model with different initial conditions 
and between different models with the same initial 
conditions. The results indicate that the range of these 
estimates is large. The mean increase over the five pro- 
jections, however, provides confidence that the number 
of RHAs will increase over Europe as the climate con- 
tinues to change. The results from the different green- 
house gas emission scenarios are somewhat inconsistent 
with what is expected based on the characteristics of the 
scenarios. The scenario A2, which is considered a high- 
emission scenario, shows lower estimates of the increase 
in RHAs than the A1B scenario, which is considered a 
middle of the road scenario. This result is from the dif- 
ferent regional climate models rather than the emission 
scenarios themselves. Up until 2050, the estimated tem- 
perature increase is actually estimated to be higher in 
A1B than A2 according to Intergovernmental Panel on 
Climate Change (IPCC). 
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To improve the estimates of the impact of heat on 
RHAs, detailed data on emergency room visits and 
admissions during summer are needed, such as the pro- 
portion of emergency department cases admitted to hos- 
pital and the death rate. A better understanding of 
vulnerable groups is needed, including how these 
groups could change over time. For example, the 
portion of the population aged over 65 will increase 
from 17% to 29% 42 by 2050, significandy increasing the 
size of this vulnerable group. COPD is mainly a disease 
of the elderly. Persons suffering from COPD appear to 
be especially vulnerable to heat, 43 and the incidence of 
COPD is likely to increase 44 with the consequence that 
the health impact of heat may increase in future years. 
This study was limited to estimating the impact on the 
total population because few age-specific relationships 
were reported and because age-stratified data were only 
available on a coarser spatial scale than the climate data, 
which would reduce the benefits of having spatial 
climate data. When spatial population data, such as the 
Eurostat data, have a finer age stratification, a more 
detailed assessment will be possible. 

Normal weather patterns in the future will be different 
than those existing today, both for average temperature 
and extreme weather events. As the number of 
heat-related RHAs is expected to increase on the sub- 
continental scale, additional national or regional projec- 
tions of the future health burden from heat are needed 
that will take into account possible changes in exposure 
and vulnerability. 



CONCLUSIONS 

Projected changes in temperature and the number of 
extreme heat events with climate change could substan- 
tially influence respiratory morbidity across Europe. 
Analyses projected that both the future proportion of 
annual RHAs attributed to heat and the relative change 
in heat-related RHAs will be largest in southern Europe, 
where they are expected to nearly triple. Eastern Europe 
can expect the smallest increase in heat-related RHAs, 
where they are estimated to approximately double. For 
all of Europe, the number of respiratory heat-related 
hospital admissions is projected to be 26 000 annually in 
the future period (2021-2050) compared with 11 000 in 
the reference period (1981-2010). The estimates pre- 
sented rely on the assumption that no additional adapta- 
tion occurs. Future studies should elaborate and 
quantify the possible effects of different adaptation 
assumptions applied to regional conditions. 
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